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•  Tortuosity  factor  of  the  3D  infiltrate  network  is  calculated. 

•  Analytical  model  for  the  tortuosity  factor  of  infiltrate  network  is  developed. 

•  Strategies  to  decrease  the  tortuosity  factor  of  infiltrate  network  are  suggested. 

•  Intrinsic  conductivity  of  the  infiltrate  network  can  be  resolved. 
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The  tortuosity  factor  is  calculated  by  solving  Laplace's  equation  of  electrostatic  potential  within  the 
numerically  constructed  three-dimensional  infiltrate  network  of  nanostructured  solid  oxide  cell  elec¬ 
trodes.  Based  on  Bruggeman  approach,  an  analytical  model  is  proposed  to  calculate  the  tortuosity  factor 
and  effective  conductivity  as  a  function  of  infiltration  loading.  The  intrinsic  conductivity  of  infiltrate 
network  can  be  resolved  from  its  effective  conductivity  using  the  analytical  model.  Good  agreement  is 
found  between  the  numerical  and  analytical  models  and  the  experimental  data  in  literature.  Parametric 
study  for  the  effects  of  backbone  microstructure,  nanoparticle  size  and  aggregation  of  infiltrate  suggests 
practical  strategies  to  decrease  the  tortuosity  factor  of  infiltrate  network. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Nanostructured  engineering  has  drawn  tremendous  attention  to 
prepare  highly  active  electrodes  for  solid  oxide  cells  (SOCs)  [1-3] 
and  gas  separation  membranes  [4  .  Infiltration  is  the  most  widely 
utilized  method  to  increase  electrode  performance  by  introducing 
nano-scale  geometric  features  where  electrochemical  reactions 
take  place.  For  specific  material  systems,  numerous  SOCs  with  the 
highest  power  density  or  the  lowest  polarization  resistance  have 
been  reported  by  using  infiltration  1,5  .  Modeling  and  experi¬ 
mental  studies  have  tried  to  link  the  performance  with  the  geo¬ 
metric  factors  of  infiltrated  electrodes.  For  instance,  the  Surface 
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Resistance  model  [6]  and  the  modified  model  [7]  calculate  the 
polarization  resistance  of  mixed  ionic  electronic  conductor  (MIEC) 
infiltrated  ionic  conductor  electrodes  as  a  function  of  MIEC  surface 
area;  The  Particle  Layer  model  links  three-phase  boundary  length 
with  the  polarization  resistance  of  the  electronic  conductor  infil¬ 
trated  ionic  conductor  electrodes  [8  .  However,  many  details  of 
infiltrated  electrodes  remain  a  mystery.  One  important  issue  is  the 
effective  conductivity  of  the  infiltrate  network  [9], 

<7eff  =  -°'0  (1) 

T 

where  gq  is  the  intrinsic  conductivity  of  infiltrate  network;  r  is  the 
tortuosity  factor;  and  0  is  infiltration  loading  (volume  fraction  of 
infiltrate  in  the  entire  electrode  including  pores).  Experimental 
observations  show  that  the  electronic  conductivity  of  electrode  can 
be  increased  significantly  by  a  small  loading  of  electronic  infiltrate 
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(e.g.  5  vol.%),  however,  can  not  be  further  promoted  to  a  high  level 
comparable  to  that  of  the  conventional  composite  electrodes 
[10,11  .  This  can  be  attributed  to  either  or  both  of  the  high  tortu¬ 
osity  factor  and  the  low  intrinsic  conductivity.  However,  the  con¬ 
tributions  of  geometric  properties  (r)  and  materials  properties  (do) 
can  not  be  well  separated  yet.  This  is  partially  because  of  the 
challenges  in  the  quantification  of  electrode  microstructures. 
Scanning  electron  microscopy  (SEM)  has  been  used  to  estimate 
infiltrated  nanoparticle  size  [7,12].  Advanced  techniques  such  as 
focused  ion  beam  -  SEM  and  X-ray  computed  tomography  have 
been  used  to  quantify  the  3D  geometric  factors  such  as  percolation 
properties  and  tortuosity  factors  for  conventional  composite  elec¬ 
trodes  [13,14  .  However,  it  is  challenging  to  capture  the  nano-scale 
features  of  infiltrated  electrodes  due  to  the  limited  resolution 
(10-50  nm).  Recently,  Kishimoto  et  al.  utilized  high-resolution 
focused  ion  beam  tomography  and  for  the  first  time  recon¬ 
structed  the  3D  microstructure  of  an  infiltrated  electrode  [15].  In 
this  work,  we  construct  the  3D  microstructures  of  infiltrated  elec¬ 
trodes  by  using  the  numerical  infiltration  methodology  of  our 
previous  work  [16  .  Then,  the  tortuosity  factor  of  the  infiltrate 
network  is  calculated  by  solving  Laplace's  equation  for  charge 
transport  within  the  3D  network.  Effects  of  backbone  microstruc¬ 
ture,  nanoparticle  size  and  aggregation  of  infiltrate  are  studied. 
Analytical  model  is  developed  for  the  tortuosity  factor  and  effective 
conductivity,  and  verified  by  both  the  numerical  and  experimental 
results. 


2.  Numerical  methodology 

The  most  common  infiltration  process  involves  depositing  a 
catalyst  precursor  solution  of  metal  nitrates  into  an  existing  porous 
backbone,  followed  by  firing  at  800-900  °C  to  form  the  catalyst 
nanoparticles.  In  our  previous  study,  a  numerical  methodology  has 
been  developed  to  simulate  the  3D  microstructures  of  infiltrated 
electrodes  [16  .  The  numerical  procedure  consists  of  two  steps.  The 
first  step  is  to  generate  backbone  matrix  by  randomly  dropping 
spherical  particles  with  an  intended  contact  angle  within  a  3D 
empty  domain,  which  is  discretized  into  voxels.  This  work  still 
employs  the  1000  nm  x  1000  nm  x  1000  nm  domain  with  a  spatial 
resolution  of  5  nm,  which  has  been  verified  to  be  sufficient  to 
represent  the  infiltrated  nano-scale  microstructures  16].  The  sec¬ 
ond  step  is  to  randomly  fit  nanoparticles  with  intended  contact 
angles  onto  the  surface  of  the  generated  backbone  matrix  one  by 
one.  The  probability  for  each  nanoparticle  coating  onto  the  back¬ 
bone  surface  is  defined  as, 

PBBS  =  wAreaBBs/[wAreaBBs  +  (1  -  w)AreaNPS]  (2) 

where  AreaBBS  [pm2  pm-3]  is  the  backbone  surface  area  after 
infiltration;  Areals  [pm2  pm-3]  is  the  nanoparticles  surface  area;  w 
is  a  weighting  parameter,  in  the  range  of  0-1.  Otherwise,  the 
nanoparticle  is  coated  onto  the  existing  nanoparticle  surface  to 
simulate  the  aggregation  of  nanoparticles. 


Fig.  1.  Visualization  of  the  infiltrate  networks  with  a  infiltrate  loading  of  10  vol.%.  Domain  side  length  is  1000  nm  with  a  resolution  of  5  nm.  Colors  from  red  to  blue  proportionally 
represent  the  electric  potential  from  1  V  to  0  V.  (a)  The  baseline  condition  with  w  =  0.9,  rNP  =  25  nm,  rBB  =  200  nm,  and  0BB  =  0.5.  b  and  c  show  the  effects  of  w,  which  is  0.5  and  0.1 
respectively;  d  and  e  show  the  effects  of  rNP,  50  nm  and  75  nm;  f  and  g  show  the  effects  of  rBB,  150  nm,  and  100  nm;  h  and  i  show  the  effects  of  </>BB ,  40  vol.%  and  60  vol.%.  (For 
interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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During  the  numerical  infiltration  procedure,  the  tortuosity  fac¬ 
tor  of  the  3D  nanoparticles  network  is  calculated  as  [9  , 


r  =  </xj0  /  (7eff 


Eq.  (3)  indicates  that  oo/oeff  should  be  a  constant.  Theoretically, 
oo/o'eff  can  be  ‘measured’  numerically  by  solving  Laplace's  equation 
for  charge  transport  within  the  3D  infiltrate  network, 


Ai//  -  0 


where  0  stands  for  the  electrostatic  potential.  To  avoid  singularity, 
the  isolate  nanoparticle  clusters  are  excluded.  Electrostatic  poten¬ 
tials  of  i//q  and  i/q  are  enforced  on  the  top  and  bottom  boundaries  of 
the  infiltrated  network,  respectively.  Then  the  conductivity  ratio  is 
estimated  as, 

o-q  _  f  \  -  V'o  A  ,  . 

<7eff  d  JJs\xpAS  ^  ‘ 

where  d  is  the  distance  between  the  top  and  the  bottom  bound¬ 
aries,  1000  nm;  A  is  the  surface  area  of  the  entire  top  boundary, 
106  nm2;  S  is  the  surface  area  of  infiltrate  on  the  top  boundary.  Eqs. 
(4)  and  (5)  are  solved  by  in-house  MATLAB  code. 


3.  Analytical  model  for  tortuosity  factor 

The  effective  conductivity  of  porous  network  is  usually  modeled 
by  the  empirical  Bruggeman  relationship  17], 


<7eff  =  <p^a0  (6) 

where  ^  is  the  Bruggeman  factor.  We  assume  the  infiltrate  network 
obeys  Eq.  (6).  It  is  noted  that  the  infiltrate  network  is  conductive 
only  when  the  volume  fraction  is  above  percolation  threshold,  </>c. 
On  the  other  hand,  the  maximum  infiltration  loading  is  the  back¬ 
bone  porosity,  </>Bb,  in  other  words,  the  maximum  conductivity  of 
the  infiltrated  network  corresponds  to  the  effective  conductivity 
when  the  pores  are  fully  infiltrated,  (7p0re-  Accordingly,  it  is  natural 
to  rescale  the  volume  fraction  in  Eq.  (6)  in  the  range  between  </>c  and 
0BB» 


^eff 


0  (pc 
_0BB  -  0c. 


4 


O'pore 


(7) 


where  the  maximum  conductivity  is  given  by, 


tfpore  —  0gg0o 


(8) 


Inserting  Eq.  (8)  into  Eq.  (7),  we  can  obtain  the  effective 
conductivity, 


°eff 


0BB 


cj)  (f)c 

0BB  —  0c. 


4 


00 


(9) 


Combining  Eq.  (3)  with  Eq.  (9),  the  tortuosity  factor  of  infiltrated 
network  can  be  estimated  as, 


T  —  (pi  0gg- 


(p  (pc 


0BB  —  0c, 

Bruggeman  factor  is  treated  as  a  fitting  parameter. 


(10) 


Fig.  2.  Tortuosity  factor  of  infiltrate  networks  as  a  function  of  infiltration  loading.  The 
effects  of  w,  rNP,  rm  and  b  are  shown  in  (a-d),  respectively.  Scatters  are  numerical 
results  and  solid  lines  are  the  fitting  of  Eq.  (10). 
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4.  Results  and  discussion 


LSM  loading,  vol.% 


C 


The  tortuosity  factor  of  infiltrated  network  is  calculated  as  a 
function  of  infiltration  loading.  Morphology  of  the  infiltrated 
network  is  governed  by  four  factors  -  the  aggregation  risk  of 
nanoparticles,  nanoparticle  radius,  backbone  particle  radius,  and 
backbone  porosity.  Thus,  the  effects  of  these  four  factors  are  stud¬ 
ied.  A  baseline  condition  is  defined  as  w  =  0.9,  tnp  =  25  nm, 
rB b  =  200  nm,  and  </>BB  =  0.5.  And  the  parametric  study  is  conducted 
by  changing  one  factor,  while  keeping  the  other  factors  unchanged. 

Fig.  1  shows  the  electrostatic  potential  distribution  within  the 
10  vol.%  infiltrated  networks  under  various  conditions.  Colors  scale 
with  electric  potential  from  1  V  (red)  to  0  V  (blue).  As  is  shown, 
when  the  weighting  parameter  decreases  from  the  baseline  0.9 
(Fig.  la)  to  0.5  (Fig.  lb)  and  0.1  (Fig.  lc),  the  infiltrated  nanoparticles 
tend  to  form  aggregates,  leading  to  insufficient  connection  between 
nanoparticles.  Comparing  with  the  baseline  condition,  the  electric 
potential  for  w  =  0.5  is  badly  distributed,  in  other  words,  the  po¬ 
tential  is  not  smoothly  evolved  from  the  bottom  side  to  the  top  side. 
According  to  Ohm's  law,  the  gradient  of  potential  is  proportional  to 
current.  Thus,  the  equipotential  clusters  in  Fig.  lb  and  c  do  not 
contribute  to  the  current  transport,  increasing  the  tortuosity.  When 
the  nanoparticle  radius  is  increased  to  50  nm  (Fig.  Id)  and  75  nm 
(Fig.  le),  the  network  is  isolated  under  the  infiltration  loading  of 
10  vol.%,  and  thus  leading  to  infinite  tortuosity.  Decreasing  the 
backbone  particle  radius  from  baseline  200  nm-150  nm  (Fig.  If) 
deteriorates  the  potential  distribution.  The  infiltrated  network  is 
isolated  by  decreasing  backbone  particle  radius  to  100  nm  (Fig.  lg). 
The  potential  is  badly  distributed  by  decreasing  backbone  porosity 
from  60  vol.%  (Fig.  li)  to  baseline  50  vol.%  (Fig.  la)  and  40  vol.% 
(Fig.  lh). 

Fig.  2  shows  the  tortuosity  factor  evolution  with  infiltration 
loading  under  various  conditions,  revealing  that  the  tortuosity 
factor  decreases  with  increase  in  infiltration  loading.  For  a  given 
loading  level,  the  tortuosity  factor  is  decreased  by  increasing 
weighting  factor  (decreasing  aggregation  risk,  Fig.  2a),  decreasing 
infiltrated  nanoparticle  size  (Fig.  2b),  and/or  increasing  backbone 
particle  size  (Fig.  2c),  consistent  with  Fig.  1.  Increasing  backbone 
porosity  (Fig.  2d)  can  also  decrease  tortuosity  factor  but  only 
slightly.  Eq.  (10)  fits  well  to  the  numerical  construction  results.  The 
Bruggeman  factors  under  various  conditions  are  in  a  range  of 
1.3— 1.7,  close  to  the  typical  value  of  1.5  for  the  oft-used  Bruggeman 
approaches  [14]. 

The  analytical  model  (Eq.  (9))  is  then  fitted  with  literature  data. 
Fig.  3  a  and  b  show  the  fitting  results  of  effective  electronic  con¬ 
ductivity  of  SrMoOs  infiltrated  network  10]  and  LSM  infiltrated 
network  11  as  a  function  of  infiltration  loading,  showing  that  Eq. 
(9)  agrees  well  with  the  literature  data.  The  Bruggeman  factors  for 
SrMoOs  infiltrate  and  LSM  infiltrate  are  fitted  as  1.3  and  1.21, 
respectively.  They  are  similar  to  that  of  the  numerical  results. 
However,  the  fitted  intrinsic  conductivities,  Go  for  the  infiltrated 
networks  are  significantly  lower  than  that  of  the  bulk  materials. 
The  intrinsic  conductivity  for  SrMo03  infiltrate  network  is  fitted  as 
170.6  S  cm'1  at  800  °C,  remarkably  lower  than  that  of  the  bulk 
SrMoC>3  material,  about  103S/cm  [10  .  In  comparison  to  the  bulk 
LSM  data  (as  the  blue  line  shown  in  Fig.  3b),  the  fitted  intrinsic 
conductivity  for  LSM  infiltrate  network  is  a  factor  of  20  times  lower. 


Fig.  3.  Evolutions  of  electronic  conductivities  of  (a)  SrMo03  infiltrated  YSZ  anode  at 
800  °C  in  97%H2  +  3%H20,  and  (b)  LSM  infiltrated  YSZ  cathode  (black)  and  LSM  -  YSZ 
composite  cathode  (blue)  at  700  °C  in  air.  Scatters  are  literature  data  [8,9  and  lines  are 
fitting  results  of  Eq.  (9).  (c)  and  (d)  show  the  calculated  tortuosity  factor  of  SrMo03  and 
LSM  networks,  respectively.  (For  interpretation  of  the  references  to  color  in  this  figure 
legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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The  low  intrinsic  conductivity  of  infiltrate  network  may  be  caused 
by  the  blocking  effects  of  grain  boundaries  and  the  impurity  of 
infiltrate  [18].  Although  it  is  difficult  to  measure  the  intrinsic  con¬ 
ductivity  of  infiltrate  network,  we  may  get  some  insights  from  the 
model.  When  the  infiltration  loading  is  near  percolation  threshold, 
the  tortuosity  factor  is  very  high,  and  thus  dominates  the  effective 
conductivity.  However,  for  typical  infiltration  loadings,  say 
10-20  vol.%,  the  tortuosity  factor  becomes  quite  small.  Fig.  3c  and 
d  shows  the  calculated  tortuosity  factor  for  the  SrMo03  infiltrate 
and  LSM  infiltrate,  which  is  3  and  1.5  at  20  vol.%  loading,  respec¬ 
tively.  In  this  case,  the  effective  electronic  conductivity  of  infiltrate 
network  may  be  dominated  by  its  intrinsic  conductivity.  On  the 
other  hand,  it  should  be  noted  that  the  model  will  always  over¬ 
estimate  the  effective  electronic  conductivity  of  infiltrate  network  if 
one  use  the  intrinsic  conductivity  of  the  bulk  material  instead  of 
that  for  infiltrate  network. 

5.  Conclusions 

The  numerical  and  analytical  models  show  that  the  tortuosity 
factor  of  infiltrated  nanoparticles  network  can  be  decreased  by 
increasing  infiltration  loading,  backbone  particle  size,  backbone 
porosity,  and/or  by  decreasing  aggregation  of  nanoparticles,  and 
infiltrated  nanoparticle  size.  The  tortuosity  factor  is  very  high 
(>102)  near  the  percolation  threshold,  while  quite  small  (<3)  at 
typical  infiltration  loadings.  By  fitting  the  model  to  literature  data, 
the  intrinsic  electronic  conductivities  of  LSM  and  SrMo03  infiltrate 
networks  are  found  to  be  significantly  lower  than  that  of  the  bulk 
materials. 
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